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I. INTRODUCTION 



Extracting a simpler representation of a dynamical system by restricting what we wish to 
describe, which is called system reduction, makes us possible to understand an apparently 
complicated dynamical phenomenon in a simple way. In particular, perturbative system 
reduction is formulated when we can express the system in question by introducing a small 
parameter. The idea of system reduction may have a long history since Boltzmann discussed 
a transport equation for systems consisting of many gas molecules. In the last three decades, 
systematic methods for perturbative system reduction have been developed and applied to 
many examples of dynamical phenomena appearing especially in fluid systems and reaction 
diffusion systems 

However, those who are not familiar with the methods often claim that they seem to 
include a sort of art. Most artistic procedure may be the introduction of multiple-scales 
variables . Actually, the choice of a set of scaled variables in the multiple-scales expansion 
requires a guess of the final result and is justified post hoc. Thus, one may wish to have a 
standard method for perturbative system reduction without any post hoc justification. Such 
a method should be completely mechanical in a sense that our intuition is excluded out. 

Here, we should remark that the multiple-scales technique was developed rather recently 
in a history of perturbation theory on differential equations. It is difficult to specify the origin 
of the perturbation theory, but Poincare was the first to discuss this subject systematically 
||. The problem was to avoid secular terms when a perturbative expansion is considered. By 
secular term is meant a term unbounded with time. Subsequently, Krynov and Bogolyubov 
proposed a practically convenient and rigorously justified mathematical method to study 
nonlinear oscillations ||. The essence of these early important works is in the notion of 
a normal form of evolution equations. That is, by making use of a change of variables, 
evolution equations can be transformed to a simpler form. This notion has no relation to 
multiple-scales. Thus, one may develop a mechanical formulation relying on the idea of 
a normal form. In fact, such perturbation theory has been discussed by Bogaevsky and 
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Povzner |7|. They developed a formalism generalizing the Poincare-Bogolyubov-Krynov- 
Mitropolsky notion of a normal form. Kuramoto also attempted to formulate a perturbative 
system reduction with a geometrical interpretation based on the similar idea ||. There 
seems to be no procedure of guessing a final form in their formulation. A subtle point is 
that the existence of a well-defined evolution equation for new variables introduced by a 
change of variables must be assumed. Precisely speaking, it may be regarded as post hoc 
justification. 

Recently, a highly mechanical formulation for asymptotic perturbation theory includ- 
ing the perturbative system reduction has been proposed by Goldenfeld, Oono and their 
collaborators PHiT|. This was called the renormalization group (RG) approach. In their 
formulation, a naive perturbative expansion is a starting point in contrast to many theories 
which reorganize a perturbation series so that secular terms do not appear. The secular 
divergence is renormalized and the RG analysis is employed. These procedures are well- 
known in field theory. The RG approach has been applied successfully to many examples. 
In the context of the perturbative system reduction, it is worthwhile noting that the ob- 
tained evolution equation through the system reduction is just the RG equation. Actually, 
recently, the Boltzmann equation |T2j and an Euclid-invariant amplitude equation have 



been derived as RG equations. One may expect that universal equations in some sense will 
be regarded as renormalization group equations. 

This paper will give a new example supporting this conjecture. We will derive, by 
employing the RG method, nonlinear phase equations describing the evolution of large scale 
modulation of spatially periodic patterns in two dimensional systems. A phase equation for 
spatially periodic patterns was first derived by Pomeau and Manneville fl4jl . Subsequently, 



phase equations for patterns have been derived in various contexts and extended to nonlinear 
ones [|TJ . Practically, a method developed by Cross and Newell seems most efficient to derive 
phase equations ]T3|. Cross and Newell also proposed a canonical form of nonlinear phase 
equations. The phase equation we will derive is equivalent to a Cross-Newell equation. 
The outline of this paper is as follows. In section [0], we introduce model equations 



to be analyzed. We analyze rather general partial differential equations including a Swift- 
Hohenberg (SH) equation |IB| . In section PIT] , a naive perturbation will be performed. Com- 
plicated calculation of secular terms will be shown in the Appendix A. In section |V| we will 
renormalize the secular terms and obtain a perturbation series uniformly valid. In section 
[V], nonlinear phase equations will be obtained by the renormalization group analysis, and a 
general formula for phase diffusion coefficients will be given. In Appendix B, phase diffusion 



coefficients for the SH eq. will be calculated. Section [Vj will be devoted to discussion. 



II. MODEL 



We consider partial differential equations (PDEs) of the form 



dw 



which describes the evolution to spatially periodic patterns in a two dimensional space, w 
denotes a scalar field, and the generalization to the multi-component cases is straightforward. 
F is a polynomial of partial derivatives of w up to the M-th order. This implies that we 
do not consider an equation with a nonlocal coupling. The assumption for F will be stated 
in order. First, F is assumed to be symmetric with respect to a parity transformation 
x — > —x, while the rotational invariance is not assumed. Second, as known in many cases, 
it is assumed that the spatially periodic solutions form a family expressed by 

w(x) = f(e,k), (2) 
e = kx + (j), (3) 

where / is a 2-7T periodic function in 9, k is an arbitrary constant vector in a certain range 
and is an arbitrary constant phase. Further, for simplicity, we assume that / can be chosen 
clS db parity symmetric one, that is, 

f(6,k) = f(-6,k). (4) 
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We can easily generalize our argument to the case that this assumption does not hold, but 
the analysis becomes complicated. 

We pay attention to a late stage in the pattern formation process, and we assume that 
almost solutions to the equation approach to a neighborhood of a family of solutions ex- 



regarded as a perturbation. Physically, this expectation is reasonable because the evolution 
of patterns becomes slower on a later stage in the pattern formation process. In order to 
formulate a perturbation problem, we rewrite Eq.([j]) as 



Do not confuse this procedure with a multiple scale expansion method. We have just set 
up the system which will be analyzed, by restricting what we wish to describe. However, 
someone may claim that the above procedure is not mechanical and too formal because e 
is not an observable parameter. We agree the former claim and will discuss it in the final 
section. Also, the latter claim may be true in the sense that we cannot say the value of e. 
We consider this equation under the asymptotic condition e — > 0+. 

The simplest example of Eq.(|5|) is a Swift-Hohenberg (SH) equation: 



where R is a control parameter. The spatially periodic solutions of the SH equation are 
expressed by Eq.(pD with Eq.(|3|) and satisfy the condition Eq.(|). (See Appendix B.) Note 
that M = 4 for the SH eq. We will develop an argument applicable to Eq.(|5|). One can 
consider the SH eq. as a concrete model. 



pressed by Eq.(^) as time goes on. Then, since the solutions satisfy F = 0, dw/dt may be 




(5) 




(6) 



III. NAIVE PERTURBATION 



We employ a naive perturbative expansion in e. The solution is expanded as 

w = Wq(x, t) + ewi(x, t) + ■ ■ ■ . 
5 



(7) 



The 0-th order solution is given by 

w Q (x,t) = f(9,k(t)) = f(k(t)x + 4>(t),k(t)), (8) 

where <fi(t) and k(t) are arbitrary functions in t. Proceeding the first order in e, we obtain 

d(f>df dk df df f , 0) 

Here, L<® is a linearized operator around the solution and formally expressed by 

(10) 
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w=f 

where 5/5w denotes the functional derivative along w. is symmetric with respect to the 
transformation x — > — x and includes partial derivative operators up to the M-th order. 

Now, we solve Eq.(H). Mathematically speaking, we should specify an appropriate func- 
tional space for the solutions, but we do not enter the mathematical details. Nevertheless, 
noting the x dependence in the left-handed side of Eq.(|9]), we can expect naively that the 
general solutions take a form 

M+l 



W l = X (l:«) a (li)> ( U ) 

n=0 

where we have introduced an abbreviation 



•E(l:n) Xi 1 Xi 2 ' ' ' X j n , (12) 



a 



(n) _ (n) 



<U- (13) 



(l:n) u- iii2---i 

We will use a similar abbreviation for tensors with indices ■ ■ -i m . (Also, dn. ) means a 
scalar variable and X(i : o) = 1.) The tensor a^! n \ is a function in 6(= k(t)x + <p(t)) and 
in t (through t dependence of d<p/dt and dk/dt) and can be assumed to be symmetric with 
respect to the permutation of indices. 

Here, in order to evaluate L^wi, we define the n-th order tensor operator l/ n ) recursively 
by the formula 

£(^-i)K#(z)] = x^L^^gix) + L ( ^ n) g(x), (14) 



where g{x) is an arbitrary function. As a concrete example, D> Q \ ZA 1 ) and for the SH 
eq. are given by 

m = R-3w 2 -(l + A) 2 , (15) 

If = -4(1 + A)d j , (16) 

Lfj = -4(1 + A)5ij - 8didj. (17) 

Then, by using the operators {L^}^^ 1 , L^Wi is written as 

M+l M+l 

#<V = E *<i*) E nC k L^ n) a^ n) . (18) 

fc=0 n=fc 

Substituting this equation to Eq.@, we obtain 

M+l 



E *(l:^(Sc) = 0, (19) 
fc=0 



where 



M+l 

Ah) _ sr^ n f(n-k) (n) fr>n \ 
°(l:fc) — 2^ "°fc-k(fc+l:n) a (l:n)' l ZU J 



n=k 



for k > 2, and 



,(1) ?(n-l) (n) 9k h df 

%:l) = E ^(2:n)^(l-n) " ( 2 1) 



1,(0) V- f(n) (») , 99 x 

(l:0) ~ h (1:n) (1:n) & ^ ■ 

Since x is arbitrary and & < - n - ) is a periodic function in x, Eq. ( |T9D leads M + 2 equations 



6jJi) = 0, (0<A;<M+1). (23) 
We can determine the tensors {a^} by solving the set of equations and obtain 

Wi = {ximnPhln + x lmPfrL + X ^)^§ + ( 3x l™Plmn + 2x ^n) j^T + ^ ■ ( 24 ) 

Here, pf^ n and do not depend on x and satisfy 



<V>!lA™ = |p (25) 

2ft<SA™ = f , (26) 

with 

A™ = \{<L^o) + (*lL$?L), (27) 

where (it, v) denotes the inner product in the space V per = {u\u G L 2 [0, 2tt), u(0) = u(2n)}, 
$o = df /dO is a null eigenvector, and $q is the adjoint null vector normalized to ($q, ^o) = 1- 
pi is an arbitrary constant and is a periodic function with a certain arbitrariness. We 
do not need to know the explicit form of a^°\ Eqs.(p4|) - (^) are derived in Appendix A. 



IV. RENORMALIZATION 

The naive perturbation series up to the first order in e is summarized as 

w = f(k(t)x + <f)(t),k(t)) + e Wl , (28) 

where w\ is given by Eq.([24]). This perturbation result does not make a sense for large |x| 
because of the singular behavior of w\. However, this result becomes meaningful locally 
around an arbitrary position X when we successfully separate the singular behavior and 
renormalize it to the arbitrariness in the zero-the order solution. This operation corresponds 
to the renormalization in field theory. 

We now show that the singular behavior in wi is renormalizable. First, we split £(i :n ) to 
X(i:n) — x (i-.n) ~ X(i:n) + X(i:n), where x^\. n ^ = X(i- n ) — A( 1:n ) is not large when we focus on the 
system around the position X. Then, the singular part in wi, denoted by wl™ 9 , becomes 

w s r = (x lmnP f2n + + x ipf ] )% + (*x*nA + 2X ^w n - (29) 

We wish to renormalize these singular terms to the arbitrariness in the 0-th order solution, 
f(kx + (f),k). Let us recall that <fi and k in / are arbitrary functions in t. We then define 
(f) R (X,t) and k R (X,t) by 
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4>{t) = R (X,t) + e50(X,t), (30) 
k(t) = k R (X, t) + e5k(X, t). (31) 

(f) R (X,t) and 8(f>(X,t) are called a renormalized phase and a counter term, respectively. k R 
and 5k are called similarly . Substituting Eqs. (|30|) and ( |5ID into /, we obtain 

f(Jsx + (f>, k) = f(k R x + (j) R , k R ) + e[(<j£x fl + + 5(j))% + Sk^l}. (32) 

afc* dk 



Comparing this with Eq. (p9|) , we obtain the cancellation condition of the singularities as 

56 + X lm J± + X lmP i 2 2 + X lP { p = 0, (33) 
5k n + 3X im p[l + 2X lP $ = 0, (34) 

where we put 

59 = 5kX + 5(f). (35) 
As a result, the renormalized perturbation series becomes 

w = f(k R (X,t)x + <f) R {X,t),k R {X,t)) 

+ 6 {[ x hnnPlmn + x hnPlm + X f (Pi ^ + ^)] 

+ (3xgA + 2xfp«)^ + aW}. (36) 

It seems strange that a singular term remains in the renormalized perturbation series. 
However, since it appears with the form x R 5ki and this term simply narrow the valid region 
of the result to \x R 5k{\ ~ 0(1). In this way, the renormalizability of the singular behavior has 
been shown, and our perturbation result is meaningful locally around the arbitrary position 
X. Thus, obviously, we can put x = X and then obtain 

w(X,t) = f(k R (X,t)X + R (X,t),k R (X,t)) + ea (0) , (37) 

where note that X can be replaced by a small letter x because X is arbitrary. This is 
a globally valid expression for w. Here, k R (X,t) and <f) R (X,t) are given by Eqs. (|30D and 
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(j3Tp with Eqs.(p3[) -(p5|). Before discussing these equations, we consider a relation between 
k R (X,t) and 6 R (X,t) which is defined by 

B R {X,t) = k R {X,t)X + <f) R {X,t). (38) 

By using 9 R , 9 = kx + <fi is written as 

= k R x R + 6 R + t(5kx R + 56). (39) 

Differentiating this equation with respect to X, we obtain 

k R = ?C, (40) 
dX 

where we have used the equalities 06/dXi = dki/dX m = and 

rf=™. (41) 
dX 

The last equality can be easily proved from Eqs.(Q) and fl3"4|). One may find that our 9 R 
corresponds to G in the Cross- Newell's formulation for the phase dynamics Actually, 
in the next section, by the renormalization group analysis, we will derive the equation for 
9 R (x,t), which corresponds to a Cross- Newll phase equation. 



V. RENORMALIZATION GROUP 

Noting k(t) = k R (0,t), we rewrite Eq . (j3ll) as 

k R (0, t) = k R {X, t) + e5k{X, t). (42) 

— * — * — * 

This equation is valid locally around X = because of the singular behavior of 8k. (See 
Eq.(|3"4"D.) Here, the capital letter X in the argument of renormalized quantities is convention. 
x can be used instead of it, of course. We now wish to obtain a globally valid expression 
for k R . The RG analysis makes us possible to do it. We first show a formal argument. 
Differentiating Eq.(fl^) with respect to X, we obtain 



10 



= -e™*, (43) 
= e(6X n p\ 3 2n + 2pf2). (44) 

Here, p^ 3 -* and p^> are functions of k(t) in their definitions Eqs.(^) and fl26|), but we are 
allowed to replace to k R (X,t) because the difference gives the 0(e 2 ) contribution in 
Eq.ffBl). Multiplying Di m {k R ) to this equation and using Eqs.([25l) and fl26|) , we obtain 

By combining Eq. (|40|) , this equation can be written as 

We expect that Eq.(f4"T|) is globally valid because it does not include a singularity. Equation 
(|4T|) is equivalent to a Cross- Newell equation |TJ|]. However, one may claim that Eq.(f4"T|) 
is valid only over the region where Eq.(T42|) is meaningful because Eq.(f4"T|) has been derived 
from Eq.(|42D. Also, the replacement from k to k R in Eq.ffHf) seems a little bit unnatural. 
We notice here that the RG analysis is not employed explicitly in the above argument. Next, 
we consider the problem from the RG point of view. 

The argument in this paragraph follows a review article by Shirkov |17|]. The RG is a 
symmetry structure with respect to the alternation of the way of giving boundary values 
|T7Hl9| • Let us see this property in our system. Q denotes a set of renormalized quantities, 
i.e. 

Q(*) = *),*«(*,*)), (48) 

where Q(X) is a function in t. From Eq.(^) and a similar expression for 9 R derived from 
Eq. (|5UD , we can write 

Q(0) = Q{X) + e6Q{X). (49) 
11 



This expression is valid locally around X = 0. Using the spatial homogeneity of the system, 
we can prove 

Q(X ) = Q(X) + e5Q(X - X ), (50) 



which is valid locally around an arbitrary position X . Then, from the form of Eq. (|50|) , we 
assume a globally valid expression 

Q(X) = F(X-X ,Q(X )) (51) 

for an arbitrary position Xq. Since Q(X) does not depend on Xo, Eq.([|l]) shows the in- 
variance property of Q(X) under the alternation of the way of giving the boundary values. 
Further, F satisfies the transitivity expressed by 

F(X - X , Q(X )) = F(X - X 1: F(Xi - X , Q(X )). (52) 

This leads to the composition law of the transformation R(X) acting on Q(Xq), which is 
defined by 

R{X)Q(X )=F(X,Q(X )). (53) 

Noting that R(0) = 1 and that R(—X) is the inverse transformation of R(X), we can 
conclude that the transformations form a Lie group. Then, according to the fundamental 
theorem of the Lie group theory pD" |, a differential equation obtained by considering an 



infinitely small transformation determines the Lie group. The differential equation is 

|S = /3(Q), (54) 
aX 



where j3 is called a generator and defined by 

dF(X,Q) 



(3(Q) 



(55) 

x=o 



dX 

Eqs.([52]) and ( |53D are called functional and differential RG equations for the variable Q 
respectively. The generator (3{Q) has been referred to as a beta function in the RG literature. 
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As seen in its definition, the generator is determined by the local property of the trans- 
formation. Thus, the locally valid expression Eq.(fl^) is enough to give the generator pertur- 
batically, and through the RG equation Eq. (|53|) , which is globally valid, we can find Q(X) 
for all X. This procedure to obtain globally improved solutions from locally valid ones are 
known as the perturbative RG method. From Eqs.(^) and (pM|), we obtain 

^ r Pf\e\k R ) = k«, (56) 

|J = k R ) = l(D-rk R )) lm ^. (57) 

These RG equations are equivalent to Eq.([|7]). The procedures Eqs. (^3|) -(^6|) should be 
regarded as a simplified argument of the perturbative RG method based on the Lie group 
theory. Recently, Kunihiro has claimed that the RG method developed by Goldenfeld, Oono 
and their collaborators is not mathematical and has proposed an envelop method against the 
RG formulation |^TJ . However, there is no necessity to consider the envelop method, though 
it may be a correct interpretation of the perturbative RG method. We recognize that the 
RG method discussed here is not purely mathematical, but still remains at a formal level 
as usual in theoretical physics. We believe however that a careful mathematical description 
will be possible. (See |22| for a mathematical formulation for phase dynamics.) 



VI. DISCUSSION 

We summarize the result: The solution is expressed in the form 

w = f(9 R ,k R ) + ea<®, (58) 

where 9 R obeys the nonlinear phase equation Eq. fl47|) , and k R is derived from Eq.(f40|). As 
far as we know, the general formula for phase diffusion coefficients has never been presented 
explicitly. When a model equation satisfies the conditions we assumed, we can immediately 
calculate the phase diffusion coefficients by using the general formula Eq. (|2"7|) . In Appendix 
B, we show the result for the SH equation. Since many variants of the SH eq. [HI satisfy 
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the conditions we assumed, this general formula is useful. However, there are some impor- 
tant classes out of our consideration. One class consists of nonlocal models which appear 
commonly in fluid dynamics because of the incompressive condition [P^j . The resultant 
phase equation is coupled to mean flow and this leads to a new type of instability called 
a skewed varicose instability JTSkPiySSfl . It seems that there is no technical difficulty to 
discuss nonlocal models, but complication will come in. The other class out of our con- 
sideration consists of models which have a different type of a family of solutions. As one 



example, since the Kuramoto-Sivashinsky (KS) equation |26] possesses a Galilei symmetry 
f27j , a family of spatially periodic solutions in the one dimensional system is expressed by 
w = f(k(x — vt) + <f),k) + v, where k, and v are arbitrary constants. Correspondingly, due 
to the existence of additional null mode, secular terms take a different form when we employ 
a naive perturbative expansion. We expect that all secular divergences are renormalized 
to k, cf) and v and the RG equation gives a correct phase equation. We remark that an 
oscillatory instability appearing in Rayleigh-Benard convection with a low Prandtle number 
is associated with the Galilei invariance [|28|| . 

A pattern breaks the spatial translational symmetry and the phase variable is interpreted 
as a Goldstone mode for the symmetry breaking. Similarly, phase equations associated with 
the temporal translation symmetry have been discussed in reaction diffusion systems |§, 
where a limit cycle solution breaks the temporal translational symmetry. The resultant phase 
equation has a similar form with a Burgurs equation. Renormalization group derivation 
of the phase equations from general reaction diffusion equations which have a limit cycle 



solution is much simpler than the derivation in this paper |29j , because the naive perturbation 



result in such systems shows that there is only one secular term proportional to t. Phase 
equations for propagating patterns will be derived similarly, though a little complication 
arises because of the mixed nature of spatial and temporal symmetry breakings. 

Let us compare the RG method with other methods for perturbative system reduction. 
The most efficient derivation of phase equations may be the Cross-Newell method. In their 
formulation, besides a standard multiple-scales analysis, a new variable is introduced 
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with playing two roles: the derivative of G gives the modulation of wavevectors, while 
has a multiple-scales relation with a phase coordinate specifying a position of the basic 
periodic pattern. That is, their formulation consists of a tactical combination of the multiple- 
scales analysis and the method of changing variables. We recognize that such an acrobatic 
procedure greatly reduces a calculation time to obtain a final form especially in complicated 
problems. Their method may be optimized for the purpose of deriving nonlinear phase 
equations. Also, the formal perturbation series can be produced systematically up to an 
arbitrary order. 

The RG derivation of phase equations need more tedious calculation than the Cross- 
Newell method. As far as we are concerned with simple problems such as nonlinear os- 
cillations, there is no significant difference between computational efficiencies by the RG 
derivation and by conventional ones. However, as shown in this paper, it seems obvious that 
a more mechanical method requires a more calculation time when the problem becomes 
complicated. Thus, one may doubt the practical efficiency of mechanical formulation. Me- 
chanical formulation seems to be apparently impractical. However, for much complicated 
problems which are intractable by manual calculation, a more mechanical method may be- 
come more useful because mechanical procedures can be programmed in computers. Further, 
mechanical methods are useful to problems for which we cannot find a suitable guess of the 
solution. 

We should remark that our derivation is not yet a completely mechanical one. As dis- 
cussed in section [H], there is subtle ambiguity where a small parameter e is put. We cannot 
formalize this procedure mechanically. This problem is not peculiar to the present problem, 
but appears in many examples, even in a quite simple example such as linear ordinary dif- 
ferential equations |jTT|. It may be honest to say that we set up systems so that the RG 
method works well. Our choice is physically reasonable, but our purpose is to eliminate such 
a phrase. In addition, let us recall that perturbative system reduction is possible only when 
we can express the system in question by introducing a small parameter. The formulation 
on this part is independent of the perturbation method and should be discussed separately. 
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However, we do not have a general way how to express what we wish to describe. Of course, 
in some cases, we can introduce suitably a small parameter with referring to the physical 
situation in consideration, but again our purpose is to eliminate such a phrase. As far as we 
consider in the present knowledge, it seems hopeless to formulate this point mechanically. 

As mentioned in section |, there are other theories aiming at a mechanical formulation. 
Bogaevsky and Povzner analyzed many examples including WKB-type problems and some 
PDE problems f7j. The correspondence with the RG method seems clear at least in simple 
problems. Their method is interpreted as systematic renormalized perturbation theory, that 
is, they consider renormalized variables from the outset. The calculation may be carried out 
mechanically, but we do not check it whether or not their formulation has the same content 
as the RG method even for complicated problems such as phase equations derived in this 
paper. 

A Kuramoto's geometrical interpretation on the system reduction comes from the notion 
of a normal form ||. In this sense, his formulation shares a common structure with Bo- 
gaevsky and Povzner. However, his main claim is that the system reduction is independent 
of choice of the method of expansion series. One may choose any expansion method as far 
as it can be performed consistently. In general, we can obtain a formal solution to a set of 
self-consistent equations under certain additional assumptions which are required from the 
outset in his formulation. According to his view, when the perturbative system reduction 
is considered, this formal solution is evaluated with referring to physical situations. This 
is one clear strategy, but seems less practical than standard conventional methods and less 
mechanical than the RG method. 

The importance of the RG method is not restricted to its mechanical nature. As is well 
known, RG has been used in diverse fields of theoretical physics in particular since Wilson 



formulated RG in a different form f50fl . Thus, by studying a common structure to several 
RG approaches, we may find a new method for system reduction. 

Shirkov reviewed that RG methods appeared in several contexts were employed based 
on the transitive property of physically relevant quantities with respect to the alternation of 
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the way of giving the boundary values |17j . He called this property functional self-similarity 



(FSS). The notion is almost equivalent to a group structure of the old style RG [p]| , |ISfl and 
is a generalization of the usual self-similarity related to power laws. In this sense, the RG 
analysis is interpreted as a method to construct a FSS symmetric solution. It is not easy to 
know the FSS symmetry beforehand. Nevertheless, if we find it before calculation, the group 
analysis [j3l] can be employed to an enlarged system consisting of the original equation in 



question and an equation describing alteration of the boundary data. Such a formulation 
was already proposed by Kovalev, Krisvenko and Pustovalov |p2|| . 

The RG method employed in this paper is identical to the old style RG developed in 
1950'. Another popular formulation of RG is the Wilson-Kadanoff scheme, in which the RG 
transformation is devised based on a physical insight of the system [|3~0| , j33|1 . Owing to its 
constructive nature, such a scheme may be useful in non-perturbative system reduction. The 



study in this direction is developing f34| . Also, related to it, the constructive RG method 
was employed to discuss asymptotic behavior of PDEs mathematically However, we do 
not yet understand a relation among these different RG approaches. 
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VII. APPENDIX A: CALCULATION OF SECULAR TERMS 



In this appendix, we solve Eq.(|ZBD with Eqs.(|2D])-(^2]). We consider the case M = 4 
which holds for the SH eq. The extension to arbitrary M(> 4) is straightforward, and the 
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result is unchanged. 

The equations in question are regarded as linear functional equations for 2ir periodic 
functions in 9 = kx + (p. In order to solve it, we consider an Hilbert space V per = {u\u 6 
L 2 [0, 27r], u(0) = u(2tt)}, where d/dxj in Z,£2 is replaced by kjd/d6 when the operator 
L^l acts on u G V per . Before solving the equations, we summarize elementary mathematical 
notions about linear algebra on the space V per . 

There exists a null eigenvector $o satisfying 

L (0) $ = 0. (59) 
In fact, differentiating the equality F\ w= f = with respect to 9, we obtain 

L {0) % = 0. (60) 

Thus, we can put 

$ = df/de. (6i) 

This null eigenvector is associated to the existence of a family of solutions parameterized 
with which comes from a spatially translational invariance of the system. 

A generalized null space in V per is spanned by only $o from the following two reasons. 
First, the other null mode which is associated with a family of solutions with parameterized 
by k is not included in the space V per . Second, if there were another null mode, a family of 
solutions would be expressed in a different form. (For a related discussion, see section |V|.) 
Let us discuss the first point. Differentiating the equality F\ w= f = with respect to k, we 
obtain 

#o)# = £<o) (f S£ + 9/ ) = 0< (62) 
dk 09 dk 

Then, it is easily found that the null mode xdf /d9 + df /dk is not a periodic function in x. 
The inner product in the space V per , which is denoted by (u, v), is defined by 



r d9u(9)v(9). (63) 
2tt Jo 



is not an Hermite operator in general. (Note that is an Hermite operator for the 
SH eq.) Then, an adjoint null eigenvector is defined by 

L<°>t${ = 0, (64) 

where $J i s assumed to be normalized as 

($S,$o) = l- (65) 

Consider a linear equation for u in the space V per 

L {0) u = b, (66) 

where b is a 2tt periodic function. Since has a null eigenvector, there is a solution only 
under the solvability condition 

«b) = 0. (67) 

When this condition is satisfied, the solution u to Eq . fl66|) is expressed by 

u = Mb + A$ , (68) 

where M is a pseudo-inverse operator of L'- -' and A is an arbitrary constant. (ML^ = 
L(°)M = 1 holds on a restricted space which consists of u satisfying ($q, u) — 0. ) 
Here, for the later convenience, we write down an identity 

~ M ^ = W/ (69) 

which is obtained by expressing Eq . (|62D in the form 



^ 1) $o + L(°)^ = 0. (70) 

Now, we derive Eqs.(|24]) - (p7|). We first show = for n > 4. We rewrite Eq.(p3[) 
with Eq.([20|) explicitly as 
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f cn'i (5) 
£ (0) "(l-5) 


= o, 


(71) 




cf(l) „(5) 
01j (5:5)"(l:5) 


+ #°>a« ) 


= 0, 


(72) 


1n f(2) (5) 
iUlj (4:5)"(l:5) 


+ 4 fW „( 4 ) 

^(4:4) "(1:4) 


+ ^ (0) "^3) 


= 0, 


(73) 


1n f(3) (5) fi *(2) (4) 
iU -^ (3:5) "(1:5) + °- tj (3:4)"(l:4) 


. qfW „(3) 
+ ^(3:3) "(1:3) 


+^ (0) "!S 2) 


= 0, 


(74) 



From Eg. (|7ip , we obtain 

aSS = pjS$ , (75) 
where is a constant whose value will be determined later. We substitute Eq.(|75f) to 



Eg . (ff2|) . Then, the solvability condition for is satisfied due to the parity symmetry, and 
we obtain a^: 

(4) (4). . - (5) df . . 

"1:4 = Pl:4^0 + 5pU^-, (76) 

where we have used Eg. (j69|) , and p^) A is a constant whose value will be determined later. 
We substitute Egs.(|75|) and ( [76] ) to Eg. (|73|) . Then, the solvability condition for a® yields 

pfl D 4:5 = 0, (77) 

where 

Ann = \^lWn^0) + W ' J^) ■ (78) 

Using the regularity of the matrix D which is assumed, we obtain 

P$ = 0. (79) 
Under the solvability condition, we obtain a^: 

„(3) _ _(3)* , zLn (4) df 

"l:3 - Pl:3^0 + 4p 1:4 ^— , {6UJ 

where we have used Eg. (|B"9"D , and is a constant whose value will be determined later. 
We substitute Egs.([75|)-(|5DD to Eg.([7|). Repeating the same argument, we obtain 
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pVI = o, 



(81) 



and 

„(2) _ Jaw , o.(3) /o 9 N 

°1:2 - Pl:2*0 + 6 P^~Q^~- ^ Z > 

One can easily see that = (n > 4) even for the case M > 4. Then, Eq.(|2"3|) with 
Eqs.(0) and becomes 

of(2) (3) , 9 f(i) _(2) , f(o)Ji) _dk h df 

(2:3) (1:3) (2:2) (1:2) "(1:1) — ^ ^ > 

*(3) (3) ?(2) (2) *(1) (1) HO) (0) _ , WW (R4) 

^(1:3)°(1:3) + ^(1:2)«(1:2) + ^(1:1)0(1:1) + ^ " ~ dt 8k t 8t 89 [ ' 

The solvability condition for yields 

6Pl^2:3 = ^ (85) 

Under the solvability condition, we obtain a^ 1 ': 

Jl) _ .(1)* , 9 _(2) /o fi N 
°1:1 - Pl:1^0 + Z P1-2qT~ ' ^° b J 

where we have used the parity symmetry and Eq.(^). p{\ is a constant whose value is not 
determined. Repeating the same argument for Eq . (|84"D , we obtain 

2p?lD, 2 = f . (87) 

is unnecessary for our argument, though we can calculate it. 

VIII. APPENDIX B:PHASE DIFFUSION COEFFICIENTS 

In this appendix, we derive concrete expressions of phase diffusion coefficients for the 
SH eq. based on the general formula Eq.([27|). A spatially periodic solution of the SH eq. is 

expressed by 

w (x) = f{6, k) = A 1 (k) cos9 + A 3 (k) cos0 + • • • , (88) 
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where 9 = kx + d>, and 



A 1 = -j=y/R-(l-k*)*. (89) 

Since A 3 ~ 0(R 3 ^ 2 ) for R — > 0, we consider only one mode. Hereafter, we will not use the 
smallness of R, but keep in mind that the validity is restricted to the small R. If someone 
wishes to obtain phase diffusion coefficients for large R, he needs to write a computer 
program. However, main features of the phase diffusion coefficients can be obtained under 
this assumption. Then, $ , df/dki, and $q are easily given by 

<5> = -A 1 sin6, (90) 
df 2(1 - k 2 )ki 



dki R-{l-k 2 ) 2 



Aicosfl, (91) 



$t = _J_ sin 0. (92) 

Substituting these expressions and {£^}^ = o given by Eqs.(|I5D and ( |T7D to the general 
formula Eq. fl2"T|) , we obtain 

n -_2n-A- 2 U +Ak k 8(! ~ k 2 ) 2 k n k m 

■L^nm Z '\- L J u nm ~ ^ n "a ryj m ^ k 2 ^) 2 ' 

Since the system has a rotational symmetry, D nm is further written as 

D nm = D,,{k 2 ) k -^ + D ± (k 2 )(5 nm - ^), (94) 

where 

l9 8(1 - A; 2 ) 2 fc 2 , s 

^ = ^- 2 - fi -(i-W - (95) 

L> ± = -2(1 - A; 2 ). (96) 
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